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LIGHT REFLECTION FROM SHOCK WAVES 


ABSTRACT 


The process by which light is reflected by a moving shock 
front is examined in some detail, and techniques are described by 
which the shock velocity, reflectivity, curvature and tilt can be 
measured simultaneously or separately without perturbing the shock 
wave in any way. The techniques developed are shown to be particularly 
applicable to the complex experimental situations encountered in 
plasma physics and give a considerable improvement over previous 
measurements in accuracy and information retrieval capability. To 
check the validity of the measurements at high pressures, a new form 
of shock tube has been developed to generate hypersonic shock waves 


in air at one atmosphere. 
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z Introduction 


1.1 Purpose and Scope of the Thesis 


The use of measurements of the scattering of a laser beam 
is now a well-established diagnostic technique in Plasma ieee. 
In the majority of cases, using quiescent plasmas, valid results 


(2,3,4) 


have been obtained, but a number of workers have reported 
results of situations where the plasma is either moving bodily, or 
expanding rapidly. In these cases pressure gradients are present in 
the plasma, and can cause refraction effects which may be confused 
with the scattering produced by the plasma. In the extreme case, 
the plasma boundary moves so rapidly that a shock wave is set up 
ahead of the plasma and a relatively strong reflection of the 
diagnostic light beam can occur at the resulting discontinuity in 
gas or plasma density. An example of this type of plasma is that 
produced by focusing a ruby laser in a high pressure gas, and 
anomalous results have been reported for the scattering of light by 


(5,6) 


such plasmas 


An investigation of the range of "scattering angles" and 
mach numbers in which light-reflection by shock waves is significant 
is thus clearly of considerable importance in the field of plasma 
diagnostics, and the exact mechanism of reflection is also of 
interest in its own right, as well as for the information on shock- 


structure which can be obtained. The primary purpose of this thesis 
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is thus to determine under what conditions light reflection can be 
expected to occur from the shock wave associated with a rapidly 
moving or expanding plasma, and how this light can be distinguished 
from that scattered by the plasma producing the shock front. Sone 
new information on the structure of shock fronts has also been 
obtained, and a powerful tool for the investigation of the general 


behaviour of shock waves has been developed. 


Le2esPreviocus Experiments 


The interaction of electromagnetic radiation with a shock 
wave has been investigated by several workers. Hey et ale 
working at microwave frequencies have used the Doppler Shift of the 

frequency of the reflected radiation to determine the shock velocity. 
The use of microwave radiation, however, is restricted to the case of 


high mach number shock waves, which produce appreciable ionization 


behind the shock front. 


An extensive investigation of the reflection of light by 


a low mach number shock wave has been carried out by Hornig and his 


(6,9. Us LL) 


co-workers over a number of years, and some useful 


information on shock structure has been obtained. A number of 


anomalies in the values of reflection coefficient derived from these, 


(12) 


and the similar experiments of Court & Mallory were observed 


(9,10) 


(11) 


however, and variously attributed to scattering of light by dust 


or water ‘ita in the shocked gas, or to shock wave curvature 
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The problemsencountered in these experiments are thus similar to 
those present in laser-plasma scattering experiments (which frequently 
involve Aen y non-planar shock reason) in that both require a 
distinction to be made between scattered and reflected light, although 
different parts of the total light-signal are required to extract 


the required information in each field. 


The experiments carried out in the course of this 
investigation follow the general lines of those of Hornig et Bis ts 
with the exception that a gas laser has been used as the light source 
instead of the intense, broadband thermal light source of Hornig et al. 
The coherence, narrow wavelength spread and low divergence inherent in 
the laser, apart from simplifying the optical arrangements necessary 

in the direct measurement of reflection coefficient, have made possible 
a number of new measurement techniques not previously available at 

such short wavelength (632.8 nm). A general outline of the 


experimental arrangement using a shock tube to produce the shock wave 


is given in fig 1.1. 
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FIG.1 


ENERAL ARRANGEMENT- SHOCK TUBE AND OPTICAL SYSTEM 
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ae Experimental Conditions 


2.1 Light Source 


The most frequently-used light source for plasma 
scattering experiments is the pulsed ruby mga but it is more 
convenient for most experiments to use a continuous-wave laser. 

In the experiments described here a Helium-Neon laser was used, 
Operating at a wavelength of 632.8 nm, and a power of 1 mW. The 
wavelength is thus very close to that of the commonly-used ruby 
laser (694.3 nm) but the power density is considerably less, and 


consequently there is much less possibility of the diagnostic laser 


beam perturbing the shock front in any way. 
2.2 Detector 


The choice of a low-power light source imposes fairly 

stringent requirements on the detector used, since the coefficients 

of reflection involved are of the order of Oss to TOee so that the 
detector must be capable of distinguishing a signal of 10a: to 10° Ws 
from the background noise, stray light etc. In addition, the response 
of the detector must be fast enough to observe the reflected tent 
signal during the necessarily short time for which the shock wave and 
light beam interact. Typical shock wave-laser beam interactions will 


have durations of only a few microseconds, so a 'bandwidth' of 


approximately d.c to 10 MHz is required. 
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The only type of light detector with sufficient sensitivity 
for use in these reflection experiments was found to be a photo- 
multiplier. The high gain R.C.A. type 8645 was chosen, and a 
preliminary investigation undertaken to determine under what operating 
conditions this type of photomultiplier could meet the requirements 


of low noise and 'bandwidth'. 


teow Eneoretical Calculation of Photomultiplier Noise and Bandwidth 


The theory of noise in photodetectors is well understood 

(13 Le ) e ° 

and need not be reproduced in detail here. The use of a 
photomultiplier entails the introduction of additional sources of 
noise, however, since this type of light-detector is basically a 
current generator (RCA 8645 tube —- maximum sensitivity, S = 1000 A/W) 
and thus requires a load resistor in order to display a time-dependent 
voltage on conventional measuring equipment. Noise is introduced both 


by the load resistor, R. and the measuring equipment itself. The type 


of measuring system used is shown in fig 2.1. 


The signal voltage generated by such a system for an 


input light power of W (watts) is then given by:- 


eee oak | (2.1) 
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FIG. 21, PHOTOMULTIPLIER CIRCUIT 
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And this signal must be observable over a total noise voltage 


comprising the resistor noise:- 


V2 = V4kr, $f (2.2) 


where 6f is the bandwidth of the system 
T, is the resistor temperature 


and k is Boltzmann's constant 


added to the measuring system noise, ves which with the equipment 
used here, (Tektronix 556 oscilloscope, type 1A5 pre-amplifier and 
P6045 field-effect transistor active probe) was constant at 1 mV 
over the entire range of signal voltages and frequency, and also 


added to the photodetector noise Me 


The Photodetector noise ies originates as a noise current 


I = VI_.e.6£.6 CSE 


P 


where G is the current gain of the photomultiplier, 


I. is the signal current 


and e the electronic charge 
clearly V_ =I (2,4) 
P pL 
= Cao) 
and To Vi/R, 


I! 


Vev_.R, .6f (2.6) 
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It is apparent from equations (2.2) and (2.6) that the 
noise voltages present are proportional to the square root of Re» 
whilst the signal voltage is directly proportional to R so that 
it should be possible to obtain any desired signal-to-noise ratio by 


simply increasing the load resistance, Re 


Unfortunately, however, the photomultiplier possesses a 
small but finite output capacitance C (in which can be included any 
stray capacitance in the tone resistor circuit, and the measuring 
probe) which couples with the load resistance to give a minimum system 


risetime:- 
T = Rc (257) 
so that the maximum frequency, f, which can be observed is:- 
Tt po (tlie oa hy Daa C2) 


and this frequency is also equal to the maximum bandwidth of the 


system, 6f, since the system must be capable of observing d.c signals. 


i Brest | 
ve SE coKRNC (2.9) 


So the quantity Rot which appears in both noise voltages is in fact 


fixed’ at 
1 . 
pad che hay ( ) 


and the value of R. used to obtain the signal voltage es (eq; 2 ci )-i8 
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determined from eq.(2.9) by the maximum bandwidth required. 


To obtain the maximum photomultiplier-measuring system 
gain and bandwidth it is thus primarily necessary to reduce the 
capacitance, C, as much as nOeetnT eS: By using connectors 
between the photomultiplier and load resistor which are as short as 
possible, and by using a low-capacitance voltage measuring probe it 
was found possible to reduce C to a total of 20 pF, comprised of 
8 pF intrinsic photo-tube capacitance, 4 pF probe capacitance and 
8 pF stray capacitance to ground. Different load resistors can then 
be employed to obtain a high sensitivity or a large bandwidth system, 
as required. The effect variation of R. on the minimum observable 
Signal power is shown in fig 2.2 which is a plot of Sf against W for 
a typical photomultiplier voltage (1200V; giving S = 600 A/W) for 
equations (2.2) and (2.6), and also for the measuring system noise 
voltage of 1 mV. Only signals which lie to the right of one of the 
solid lines on the graph can be distinguished from the appropriate 
noise signal, (signal to noise ratio of one). In practice, of course, 
slightly smaller signals can be observed, as shown by the experimental 
points on the graph, which were obtained by using neutral density 
filters to vary the laser input power to the phototube with a number 
of different load resistors, Re The minimum observable power is 
found to correspond to a ratio of signal to total noise of about one 


to ten. 
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FIG. 22 BANDWIDTH VS. LIGHT SIGNAL 
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The graph shows that the system is capable of reducing 
the resistor and photomultiplier noise components to negligible levels, 
leaving the noise in the measuring system. Even with this relatively 
high noise level, however, system bandwidths of 10 MHz are attainable 
even for input powers as low as ‘ie cad heme watts, so that using such 


a measuring system, reflection coefficients of tor to ot should 


be measureable, using a 1 mW laser as the light source. 
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Ler Theory 


3.1 Summary of Shock Tube Theory 


Syed latcroduction 


The production of shock waves, and the operation of the 
conventional shock tube is now well-understood, and is discussed in 
some detail by ieee Only a brief outline of the theory will 


be given here. 
3.1.2 Compressible gas flow 


Three equations can be conveniently set up to describe 
the motion of a compressible gas in one dimension relative to the 


laboratory frame of reference. 


a) The equation of conservation of mass 


3p 1, Fo. a 
A ye + aes (puA) 0 


where 0 is the gas density 
u is the gas velocity 


and A is the cross-sectional area of the mass of gas considered. 


b) The equation of conservation of momentum 


i 
ot il ox De dat 


where op is the gas pressure 
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cy The equation of state of the gas 


p= Roet PA) 


where T is the gas temperature 
and Ro is the universal gas constant 


Assuming that we have an ideal gas. 


If we further assume that the gas is contained in a tube 
of constant cross-section, A, and that the propagation of the pressure- 
Beeturpancesis adiabatic, equations (3.1), (3.2) and (3.3), can be 


combined with the relationship 


pp’ = “constant” (3.4) 


to give the propagation velocity of a pressure-pulse through a gas 


as 
r y-L 
: Vee Dae, 2 
=¢, (e - 2 3.5 
4D fe) APS om Y-1 ( ) 
where C, is the velocity of sound 
and Py the pressure in the undisturbed gas 


3.1.3 Formation of Shock Waves 


It follows from eq.(3.5) that for a weak pressure pulse 
(p=p_), U_ tends towards C_, the sound speed in the undisturbed gas, 
O Pp fo) 
but that for large overpressures, p, the disturbance may be propagated 


considerably faster than the velocity of sound. Furthermore, it is 
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evident from eq.(3.5) that the Ol increases as the pressure, 

p, increases, so that the higher pressure region of any pulse will tend 
to propagate faster than the preceding lower pressure regions, causing 

the pulse to 'steepen up' as it propagates, and ultimately to form a 


pressure discontinuity - termed a 'shock wave'. 
.3.1.4 Shock Relations for an Ideal Gas 


The equations derived in section 3.1.2 for a simple 
pressure pulse no longer hold once the shock wave has developed, 
since the assumption of adiabatic (more exactly - isentropic) 


behaviour at the pressure discontinuity is no longer valid. 


We consider a frame of reference moving at the shock 
velocity Vs" Distinguishing quantities in the undisturbed and 
shocked gases by the subscripts 0 and 1 respectively, we define for 


convenience, three new quantities:- 


Shock Mach Number M = vein (3.6) 
Shock Strength y= P,/P, C3 ie) 
Shock Compression a = p,/e, (3.8) 


The law of Conservation of Matter then gives:- 


“= es = 3.9 
pov, = 0, (¥,-¥) =m (3.9) 


where m= the mass of gas "swept-up" by the shock front in unit time 


and u, = gas flow velocity behind the shock 
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Newton's 3rd Law of Motion gives:- 


The Law of Cons 


ervation of Energy gives:- 


gf 2 a Je 2 
Pe. + mU 5 + amv. = Py (Vv, u,) + mU, + am(v uy) 


2 


where U is the internal energy of the gas. 


(3.10) 


Gay) 


Eqs. (3.9), (3.10) and (3.11) are known as the Rankine-Hugoniot Equations 


For an ideal ga 
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so that U = 
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Number and Shock Compression as a function of shock 
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3.1.5 Shock Tube Theory 


The typical shock tube consists of two sections, separated 
initially by a gas-tight diaphragm -a driver section, at a high 


pressure, P,, and an expansion chamber, at a lower pressure, Po: 


When the diaphragm is ruptured, a shock wave is formed, 
and propagates down the expansion chamber, whilst an expansion wave 
propagates back through the driversection. In order to avoid 
discontinuities forming in gas density and pressure behind the shock 
front, it is clear that the gas flow velocities behind the shock 
and the rarefaction waves must be identical. Since it is possible to 
have different gases, initially in the drive and expansion sections, 
we denote quantities referring to the region behind the shock front, as 
far as the 'contact surface' which separates the two gases, by the 
subscript 1, and quantities referring to the drivergas between the 
contact surface and the rarefaction wave by the subscript 2. The 
subscripts 0 and 3 refer to the undisturbed expansion chamber and driver 


section gases respectively. 


From eqs.(3.1), (3.2) and (3.3) we can derive the gas 


flow velocity behind the rarefaction wave as 


u. = -—— C215 —1}p | (3.20) 
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whilst the Rankine-Hugoniot equations can be solved in terms of 
the gas flow velocity behind the shock wave, to give 


oe C (1-1) (Py /P 7D 
(2 aes (3.21) 


—— A ee 


ov Berke 
\ (1+u 5) (py /P, + uy) 


And, to avoid discontinuities at the contact surface, we must have:- 


Bana ee (3.72) 
us =u, (3.23) 
and we define y = P,/P (=p,/P3) (3-24) 
so we have:- 
Y371 ; 
Sy IOUNAY is aiemarsa A eo aes Ale (3.25) 
C,y Ui) Grn) y371 | P3 . 
Lo. 


which defines the shock strength, y, in terms of the initial pressure 
ratio Pi /P3: and the quantities Uo» Y3? Co» Cy, which are characteristic 


of the gases used. 


We can thus determine the shock strength, y, for any 
given experimental conditions, and then use eq. (3.17) to predict: the 
Mach number of the shock wave produced. Fig. J eleiewastyo.cal splot 


of y vs P3/Po3 whilst fig. 3.2 shows M as a function ap he 
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FIG. 31 SHOCK STRENGTH VS. DRIVING PRESSURE RATIO 
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FIG. 3-2 MACH NUMBER_VS, SHOCK STRENGTH 
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3.1.6 Production of Strong Shock Waves 


It is apparent from eq.(3.17) that a high mach number 
shock requires a large value of y, the shock strength. Consideration 
of eq.(3.25) shows that the largest values of y are obtained if the 
pressure ratio P3/Po is made very large. The right-hand side of 


eq.(3.25) then simplifies giving:- 


ue (1-9) Gy-1) Ue G46) 
GVcn Gt) oP 
or uf 2y Yt) (Ay eee 
(Yg-1)° oF 
ify >> fi. 


It is thus essential, once P3/Po has been made as large 
as possible, to choose the gases so as to give as large a sound-speed 
ratio (C,/Co) as possible. For this reason, Helium is commonly used 
as the driver gas, and air as the expansion chamber gas (C,/Cy = 2.91). 
Fig. 3.1 shows the variation of y with initial pressure ratio P3/Po 
for a Helium-Air shock tube. Fig. 3.2 is a plot of eq.(3.17) for air, 
whilst fig. 3.3 is derived from 3.1 and 3.2, and shows the Shock Mach 
number as a function of P3/Po: Teas evident trom fig. 3.Jethat for 
shock waves moving into air at 1 atmosphere, the driver section pressure 
becomes prohibitively high for shock waves beyond about Mach 3.5, and 
the simple pressure-driven shock tube described here cannot be used 


to produce strong shock waves in atmosphere pressure air. 
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3.1.7 Strong Shock Equations 


In addition to the simple relation for y (eq.3.27) the 


following approximation can be made, for the case where y is large 


(strong shock waves). 
Ju a ae 
(from eq.3.20) 


and from eq. (3.17) 


cat CP kanes 
~ \ Lu, 


giving M = \/(y +1) a 
: re) 
so using eq.(3.26) we find 


Yt1 ¢, 
Ne A Coamee 4) oe) 
is 5 


for strong shock waves. 
i.e. once the initial pressure ratio P3/Po is large 


a strong shock wave (y>>1), the mach number becomes 


(3.28) 


(3529) 


(3.30) 


(36.31) 


enough to produce 


independent of the 


pressure ratio, and can only be further increased by increasing the 


sound speed ratio, C,/Co. 


It should be noted here that the value of M given by 


eqns. aL) is an absolute maximum, and can never be achieved in practice. 


Eq. (3.31) would give M = 10.5 for Helium-Air for instance, at which 


value the shocked gas would be heated sufficiently so as to become 


ionized, so that Y_ would no longer be a constant, and the equations 
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must thus be modified to account for the non-ideal gas behaviour of 

a strongly shocked gas. In addition the shock wave will only achieve 
its full theoretical strength if the driver gas can escape from the 
driver section fast enough. The maximum speed of the driver gas is, 
of course Uy» and is given by eq.(3.42) for the strong shock case. 

The value of a high C, is again immediately apparent. This effect 

can be mitigated somewhat by the use of a non-uniform driver tube 
cross section (converging conically towards the diaphragm) to increase 
the gas efflux velocity. An increase of 10% in M has been obtained by 


(16) 


Alpher & White by this method. 


High values of C, are attained in practice by heating the 


(17) 


3 
driver gas, either by combustion 


(18,19) 


or, more recently, by means of an 
electrical discharge In chapter 7 a form of shock tube which 
combines a coned driver section, helium filled, with intense heating by 


means of an electrical discharge to achieve Mach numbers greater than 


ten into air at one atmosphere is described. 


3:2. Theory of Light Reflection from a density discontinuity 


Bezaie First order theory 


(7) 


Cowan & Hornig have shown that to first order, the 


reflectivity of a shock front is given by:- 
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where n(x) is the refractive index profile through the 
shock front, x being the spatial coordinate 
parallel to the axis of the shock tube. 


X is the wavelength of light used. 


) =n/,-i is the grazing angle of incidence. 


v-1 


Wl 


and aI 


Equation (3.31) has been found to give values for R 


(10) 


larger than are found in practice » and it is apparent that the 


simple first-order theory used is not adequate, at grazing incidence. 
So2.2.) pecond order Effects 


At grazing incidence refraction, attenuation and multiple 


reflection effects may become important. Pert & Shige 


have analyzed 
these effects in some detail, but for these experiments a more 


approximate analysis proves to be sufficient. 


It is clear that the second-order effects all tend to 
reduce the value of the coefficient of reflection. For the case where 
the optical path length of the reflected beam is small compared to the 
wavelength of the radiation used, the shock front will behave as a 
planar refractive index discontinuity. From classical optical aWeOry Sane 
the intensity reflection coefficient is given by the Fresnel relations, 
which reduce, at grazing incidence to:- 
aaeP (3.33) 
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where 6, = wh- | (37734) 


Oey 

fe) b 5 (355) 
being the grazing angles of refraction and incidence, respectively. 
Eq. (3.33) which ignores second-order effects thus defines an upper 


limit,R,, to the true value of R, 


In general, however, the optical path length, A, is 
not sufficiently small to give zero phase difference between 


Successive multiply-reflected beams, i.e, 


r 
A = 2d cos 85 rs (3.36) 
where d = incident beam width 


The multiply reflected beams can thus destructively 
interfere, as the angle ae is reduced, (at larger path lengths) 
and in the limit, complete interference occurs, and the contribution 
of the multiply reflected beams can be neglected, in determining a 
lower limit to the reflection coefficient. Absorption of light is 
also negligible in most experimental situations, so the only 
contribution to be included is that due to refraction in the shock 


fron. 


The effect of refraction may be determined by summing 


the contributions to the electric field vector, E, from elements of 
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thickness, 6x, of the shock front. The decrease in the electric 


vector is thus given by:- 


EE = BOE, = ESE (Soa) 


since the phase difference between the elemental components will be 


small, at grazing incidence, and the fields can be added scalarly. 


Applying the Fresnel relation (3.33) to the element, 6x, 


of the shock front then yields a deflection in beam direction given 


Dye — 
33 ea (3.38) 
eee 20 ; 
E 6 
Pee a ak 
so that ln eae a9 ln F (3.39) 
fe) ) 
since the incident wave has electric vector strength ES and angle BA 
and the transmitted wave has strength Ej» at angle 8° BaseG3.o7) 
and (3.39) then give the reflected field intensity as:- 
er 
ad iC La as 
E= £5 1/2 OE 
IS) 
1 
so the intensity reflection coefficient will now be:- 
Key 
(apf ?-oid! ) 
Ras (3/41) 
2 8, 


defining a lower limit to the value of R. 


R, <R< Ry (3.42) 
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(20) 


The more careful and detailed analysis of Pert & Smy gives an 


identical expression for R,, and also a general relationship for R, 


Diaz 


which is complex and cannot be evaluated analytically. Fig. 3.4 
shows, however, that R) and Ry do not differ significantly in these 
experiments, and either relationship (3.33) or (3.41) can be used 
in practice. 


The angle 6, in (3.33) and (3.41) is determined from Snell's 


th 
law: 
n cos 8 = (ntén)cos 84 (3743) 
where 6n is the refractive index change across the shock front. For 
small angles, 8? 81> (3.43) can be expressed as 
of a be aoe | (3.44) 
1 fe) n 
and=substitution of (3.44) in (3733) ox, (3°41) %gives, 
ais 
Rese OY? Ap (3.45) 
46 s) 
o Oo 
which reduces to: 
Re CB? (3.46) 
46 
) 
(8) 


as found by Cowan and Hornig and shown in Fig. 3.4 ('first order 


theory'). 
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; : 4 ; ia 
It is thus seen that R varies as 1/6, » near grazing incidence, 
and any change in inclination of the shock front to the tube axis, due 
to curvature or tilt of the front will have a large effect on the 


measurement of R. 
Bazeo Calculation mr Reflection Coefficient R 


Any quantitative calculation of R requires a knowledge of 
the refractive index change, $n, across the shock front. For a gas, 


at any given wavelength, the refractive index is given bee 


n = 1+Kp (SEs) 


where K is a constant, known as the Gladstone-Dale constant. The 
; rf (15) 
value of K for air, at visible:wavelengths is 0.226 . 


For the gas in front of, and behind the shock front we then 


have:- 


n 


1 + Kp (3.48) 
Oo Oo 


n 1 + Kp, (3.49) 


1 
so that n.-n = 6n= K(p, 7°) 


Il 


or én Kp ,(n-1) (3750) 


since n the shock compression is defined by 


p 
in = _i (350) 
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n is given by eq.(3.18):- 
d Bary, 
1” THiy 
ee 
where Hho ona ae ((3.19)) 


and y is related to the Mach number by eq.(3.17) which gives: 


2 
y = (lt) = 


Somtnat by using eqs.(3.52), (3.53) and (3.54) the quantity 


22 can be determined in terms of the experimental variable (M), 
O 
and the constants of the shocked gas (n)3 uy)» giving:- 
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From eq.(3.43) we have 
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so that 81 can now be determined (6, is an experimental variable), 


and hence Ri or R, can be calculated from eqs. (3.33) and (3.41) 


respectively. 


The first and second order approximations 


[eqs. (3.46) and (3.45)] can be calculated directly from the value 


of on 
n 
O 


The 
angle is shown 
experimentally 


approximations 


theoretical .varfation otf Ry with mach number and 


1s. 3.) eande3;6, | [he valuerot R, is 


indistinguishable from R,, (fig 3.4) and the 


dys 


to R are of little interest in themselves, and are 


therefore not shown. 
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FIG. 35 REFLECTION COEFFICIENT VS. MACH NUMBER CTHEORY) 
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FIG. 3:6 REFLECTION COEFFICIENT VS. GRAZING ANGLE (THEORY) 


FOR A SHOCK WAVE MOVING INTO AIR AT 1 ATM. 
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IV Simple Measurement of Shock Velocity and Reflection Coefficient 


4;). introduction 


This initial experiment followed directly the method of 


Cowan & Tay except that the laser-photomultiplier system 
described in chapter two was used in place of the white-light 
system and complicated optical arrangement of previous workers. 


The experimental arrangement is shown in general outline in 


figure 1.1 and in more detail in figure 4.1:- 
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4.2 Shock Tube 


The shock tube used in these experiments was a low 
mach-number pressure-driven tube 51 mm in internal diameter and 
6000 mm long. The shock wave was generated by rupturing a 0.13 mm 
thick 'mylar' diaphragm by overpressuring a 930 mm long ‘driver section' 
at the upstream end of the tube. Typical rupturing pressures were 
2x10° yn (5 atmospheres). The theory of operation of shock tubes 
of this type is well-understood (see chapter 3). In these 
experiments Helium was used as the ‘driver’ gas, with air at one 
atmosphere as the ‘shocked’ gas. Under these conditions shock 
waves having Mach numbers of up to approximately 3.5 in be obtained 
(chapter 3). Approximately 1000 mm of tube length were allowed 
downstream of the diaphragm for the shock wave to develop and 
stabilize, and all measurements were made in the following 1000 mm 
of the tube, the overall distance of 1700 mm from the diaphragm 
being sufficiently short at Mach numbers above about 1.4 to avoid 
interference between the shock wave and the rarefaction which 
propagates back from the ruptured diaphragm, to reflect off the 
back of the driver section, and thence down the shock tube, at 
approximately the speed of sound in the driver gas. "Mylar' was 
used as a diaphragm material since its mode of rupture is a 
simple arid tine across a single diameter, so that no fragments of 


diaphragm material are introduced into the shock tube. 


The velocity of the shock wave was measured by observing the 
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output from a number of piezo-electric transducers (‘pressure probes') 
mounted flush with the shock tube wall at various points along the 
length of the tube, on either side of the test region. The pressure 
probes had a risetime of about 10 usec and an output of approximately 
one volt under the experimental conditions giving a sufficiently 
accurate velocity determination over separation of 1000 mm as used. 
The average velocity derived from 14 measurements taken over 1000 mm 
lengths of the shock tube on either side of the test region, 
simultaneously with the reflection measurements was (680 + 80) m/sec; 


or Mach (1.8 + 0.2) at room temperature. 


4.3 Optical System 


Light was introduced to the shock tube, and the resulting 
reflection observed by means of a window section having two 200 mm 
by 10 mm flat plate glass windows 10 mm thick, mounted in the 
shock tube wall parallel to the axis of the tube, and diametrically 
opposite each other. The windows were recessed by 2 mm at their centre- 
line, and 5 mm of the edges, so care was taken always to work as close 
as possible to the upstream end of the windows so as to minimize the 
effect of any disturbance to the shock wave introduced on traversing 
the wall-discontinuity introduced by the windows. The laser beam 
passed through entrance and exit slits at the windows to define the 
angle of incidence, i, of the beam as shown in figure 4.1. The method 
used to calculate the shock velocity, and the effect upon the 


measurements of any tilt or curvature of the shock wave are shown in 
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FIG. 42 SHOCK VELOCITY MEASURMENT FROM PULSE DURATION, 
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FIG.43 EFFECT OF SHOCK TILT 
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figure 4.2 
4.3.1 Method of Measuring Shock Velocity from Duration of Reflected Pulse 


The ‘interaction distance’, 1, for which a reflected pulse 


can be observed is given by: 


tee 2iaCObe 2 (fan) 


so the duration ty of the reflected pulse is 


th = RE eot 1 (4.2) 
where tenes shock tube radius 
a shock velocity 
i = angle of incidence 


4.3.2 Effect of Shock Tilt (fig 4.3) 


For a shock tilt of ¢ (fig 4.3) the interaction distance l 


is reduced by an amount jae tan ¢, giving 


1 


2r, (cot i - tan 6) (4,5) 


i 


re 
and t 2 — (cot i - tan $) (4.4) 
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pie SEE OF SHOCK CURVALURE 
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Ao oeeeettect of Shock Curvature 


The effect here is not to shorten or lengthen the duration 
of the reflected pulse, provided that all the reflected light is 
collected by the photomultiplier, but to alter the angle through which 
the incident beam is deviated by an extreme amount of 2 nts 
where ts is the radius of curvature of the shock front. The ee 


used here is not sensitive to such changes in angle. 


The shock velocity determined from the pulse length, averaged 
over 10 results at each of seven different angles of incidence was 
found to be (800 + 200) m/sec, independent of the angle of incidence - 
in good agreement with the value obtained from the pressure probes, 
but of considerably less accuracy. The relatively large variation in 
the velocities determined from the pulse lengths was probably due to 


a@erandom Variation in shock: tilt. 


Theory (chapter 3) indicates that the reflection coefficient 
should vary approximately as the fourth power of the grazing angle of 
incidence, we Fig 4.5 shows the measured intensity reflection 
coefficient as a function of angle of incidence, together with the 


theoretical curves derived from eqs. (3.33) and (3.41). 


Although the experimental error was fairly large for these 


preliminary measurements, fairly good agreement with theory is 
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obtained, particularly at the smaller angles. The results thus indicate 
that the general procedure used gives valid results, and that the measuring 
system performs sufficiently well. A more accurate and sophisticated 
investigation of shock reflectivity was therefore undertaken, following 


the same general lines. 
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RIG REFLECTION COEFFICIENT VS. ANGLE CEXPERIMENT) 
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V Measurement of Shock Velocity, Curvature & Reflectivity using a 


Differential Interferometer 
ee eae SN LeLLeLOmecer 


5.1 Magnitude of Doppler Shift 
(7) 


Microwave measurements have established the validity 
of using the Doppler Shift in the frequency of electromagnetic waves 
reflected from a shock front to determine the velocity of the shock 
wave. The use of long wavelength radiation gives poor spatial 
resolution, however, and optical wavelengths appear very attractive 
from this point of view. Their use has only become practicable in 


recent years, however, with the advent of the gas laser, with its narrow 


bandwidth and long coherence length. 


The Doppler Shift, Av, is given by:- 


ee cos i (5.1) 
v c 
or Av = & cos i (5:2) 


where v is the frequency 


and X is the wavelength of the electromagnetic wave 


Thus for a shock wave at Mach 2 (in air) using a HN, laser, 


2 x 664 


5 cos i (5.3) 
633 x 10 


Av = 
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a fejalep 2d OnceacliOe, 


cos i (5.4) 
so that even for grazing angles of ia or so, frequency shifts of 
the region of 1 MHz are expected. However,earlier results (chapter 
4) Andicatedtithatsar shocks front) tilt of 1° or more might be expected 
and that difficulties with stray light reaching the photodetector 
are experienced at angles less than about 2°. These considerations 


limited the maximum grazing angle of incidence to 3.5 x re radians 


giving a minimum Doppler shittiofs 34, Muz, 
5.2 Observation of Doppler Shift 


Any small frequency shift in a beam of electromagnetic waves 
relative to some beam of standard frequency is most conveniently measured 
by observing the heterodyne frequency produced when the frequency- 
shifted beam is mixed with a part of the reference beam, the 
heterodyne frequency being equal to the frequency difference 
between the two beams. Only part of the reference beam can be used, since 
it is necessary to operate the photomultiplier near saturation. 

This technique suffers from two major drawbacks, under the 


present circumstances :- 


(a) The incident, (reference) beam is three to six orders 
of magnitude more intense than the beam reflected off the shock front 
so to get efficient mixing to produce a clearly defined heterodyne 
frequency, for which beams of similar order of intensity, 
are required, a variable, but extremely dense calibrated filter 


would be required. Such filters are unobtainable commercially and 
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and can only be roughly approximated to by the use of absorbing 
liquid solution (e.g. CuSO, ) of variable concentrations. 

(b) The reflected beam is not fixed in point of origin 
so that the beam moves bodily to one side as the shock wave 
traverses the interaction region. This lateral movement can be 
transformed into an angular one, by means of a suitable lens 
system, but the process of mixing the incident and reflected 
beam to produce the heterodyne signal then occurs with a 
continuously varying angle between the beams, with a resulting 
variation in beat frequency. To overcome both these problems, 
the Differential Doppler Interferometer described below was 


developed. 
5.3 The Differential Doppler Interferometer 


The general arrangement of the interferometer is shown 
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GT DIFFERENTIAL DOPPLER INTERFEROMETER 
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The optical arrangement is symmetrical about the 
centre-line of the shock tube. The light beam from the laser 
is first directed onto a beam splitter, having a 50% reflecting 
front surface, and a 100% reflecting back surface, and a thickness 
of approximately 8 mm. Two beams of aeveostins eal equal intensity 
are reflected from the splitter, one for each surface, and by 
varying the angle of the splitter relative to the incident beam, 
the lateral separation of these two beams can be varied. The 
two beams are brought to a common focus ve the point within the 
shock tube where the shock velocity is to be measured, using a 
lens of about 200 mm focal length. As the angle of the splitter 
relative to the beam from the laser is changed, the angle, 66, 
between the two beams incident on the shock front alters, typical 
values being in the range 0.5 - 1.5 degrees (for 483 a 10° 
‘radians The common angle of incidence, 2° (or beam-shock wave 
angle 8) may also be altered, by adjusting the position of the 
splitters relative re the lenses. To obtain easily measureable 
reflected signals angles of incidence in the range 80°-87° were 
used. 

The reflected beams from the shock wave are brought 
parallel again by means of an identical lens to that on the input 
side, and recombined at the front surface of a second "beam splitter". 
Both beams are Doppler-shifted in frequency, on account of their 
reflection from the moving shock front, but by differing amounts, 


due to the difference, 668, between the angles of incidence. A 
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small frequency difference thus exists between the two reflected 
beams, and is observed as a heterodyne or 'beat' frequency in the 
recombined light reflected off the second splitter. 

Because of the very small inter-beam angle, this 
frequency can be made very much less than the conventional 
Doppler frequency (section 5.1). In addition, the intensities of 
the two reflected beams are very nearly equal, since the angles 
of incidence are only very slightly different, so that. a large 
depth of modulation of the output beam can be obtained, and both 
reflected beams undergo the same lateral movement as the shock 
wave traverses the interaction region. Both the major difficulties 
inherent in the conventional Doppler shift measurement are thus 
eliminated by the use of the Differential Interferometer. It is 
also apparent that any tilt or curvature of the shock front, up to 
angles comparable with the grazing angle, 6, will still result in 
an output being observed, albeit from a slightly different point 


Oneness Axis. OL ther shock, tube, 


5.4 Modulation Depth 


In general, three distinguishable beams are reflected 
from the second splitter, for front surface reflectivities of 
40 to 60%, as shown in fig. 5.1 at A,B and C. It is immediately 
apparent that beam A has no component originating from reflected 
beam (2) (that furthest from the downstream shock tube axis) and 


will thus be unmodulated. Beam C is produced by multiple 
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reflections at the second splitter, and is thus of considerably 
lower intensity than the beams A and B. Further beams, produced 
by still higher-order reflections also exist, but, are generally 
of negligible amplitude. 

An elementary consideration of the means by which the 
two splitters produce the various beams A,B,C....... » shows that 
a common splitter front surface intensity reflectivity of 38% 
will produce equal intensity components in beam B, from the two 
reflected beams 1 and 2, resulting in 100% modulation of this beam. 
Under these conditions the relative intensities of the beams 
A,B, and C will be 3:6:4, with the higher-order beam intensities 
negligible; and some 30% modulation can be expected in beam C. 
To obtain maximum depth of modulation in the reflected signal, a 
stop may be inserted in front of the photodetector, so as to 
admit beam B only. This epePthte limitation severely restricts 
the time for which the reflected signal is received, however, and in 
practice all the light from the second splitter was allowed to fall 
on the photodetector surface. It can easily be shown that a 
maximum modulation of 65% (which should be readily detectable), can 
be expected, under these conditions, Apne Se 35% of the light 
from the laser being lost in the two splitters, for a reflectivity 


Of 387. 


5.5 Measurement of Shock Velocity 


Equation (5.2) gives the frequency shift in a single beam 
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reflected at an angle of incidence, i , from a moving shock front, 


of velocity v5 2s 
Oe Cee ’ 
vi = y cos i (555)) 


So the differential beat frequency generated by mixing two beams 


differing in angle of incidence by an amount di will be 


he ete col nana (5.6) 
Vv 
=D mame os i - cos(itéi) ] C5 ah) 
‘s 6i 
= 4 ; Sin sh we 7 copes y 
Ys 
or Vres 2 ne é6i sini (5.9) 


assuming 6i<<i, and is so small that sin di = di 


Because of noise considerations, adequate signal strengths for 
measurement are only obtained at very large (=1/,) angles of 


incidence, so it can also be assumed that sin i = 1 giving 
Sie. 
v. = 2 aos (530) 


so the shock velocity can be readily obtained, since is known, 
and 6i is fixed by the geometry of the system. Equation (5.10) 


may also be written as 
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Nes ee éi | Gealay: 


where N (= t Vy) is the total number of beats observed and 


D 
1 & Ve:ty) is the length of the region in which the shock 


wave - light beam interaction occurs. 


The variation (if any) of shock velocity across the shock front 
during its transit of the interaction region can clearly be 
determined by measuring the duration of each cycle of the observed 
waveform, since successive time-segments of the waveform represent 
reflection from successive lengths (of \ 3 mm) of the diameter of 


the shock front which lies in the plane of the optical system. 
5.6 Measurement of Shock front curvature 


The initial alignment of the differential interferometer 
is such that the two beams are focused on the axis of the shock 
tube, so that in general, as the shock front traverses the 
interaction region, the two incident beams are not reflected from 
the same point in the shock front. Any curvature of the shock 
front will then mean that the angles at which the reflected beams 
emerge will differ from their respective angles of incidence, 
except when the reflection occurs from a point on the axis of the 


shock tube. The two beams will thus be displaced relative to each 
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other at the second splitter, and there will be a degradation of 
the mixing process, with a consequent reduction in the depth of 
modulation of the signal detected by the photomultiplier. 

| From equation (5.11) and the geometry of the optical 
system it can be shown that the shock front radius of curvature, 


r_, in the plane of the interferometer is given by:- 
rs 


1e fe tapes: — (o9t2) 


2(8+56)6 . 
where N is the number of beats observed 


dis the diameter of the unfocused laser beam 


In practice, d is not a well defined quantity, since the laser 
beam has a circular cross-section, not allowed for in the two- 
dimensional analysis used in deriving (5.12), and furthermore, 
has a non-uniform intensity distribution across that cross- 
section. An effective value of d can be obtained, however, by 
calibrating the system with a reflector of known radius of 
curvature, and thus curvature of the shock front can be 
calculated from a count of the number of beats observed in the 
output pulse. If beats are observed for the entire duration of 
the observed pulse, only a lower limit on the value of the 


radius of curvature can be determined, since there is always the 
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possibility that a system with a longer measuring time would 


have given a finite number of beats. 
5./ Measurement of Shock Reflectivity 
Siu lin Rejection. of. Stray and Scattered Light 

If the interferometer is set up as described in sections 
5.3 and 5.4, a number of 100%-modulated beats will be observed, 
during some part of the output pulse, provided that the only 
light entering the detector is the two beams reflected from the 
shock front. In practice, however, there is always some stray 
room light present, in addition to the possibility of light 
being scattered out of the main beams by dust or water vapour 


(910,12) or from the shock tube 


particles in the unshocked gas, 
walls and window sections. This light, however, has no fixed 
phase relation to the input beams, and hence to either output 
beam, and therefore does not give rise to interference at the 
second splitter, but is collected by the photodetector, and 
observed as ad.c. level, with the beats due to the true 
reflected beams superimposed. 

The amplitude of the beats observed in the total output 
is thus equal to the combined intensities of the two reflected 
beams, so that if the total intensity of the input beams is 


measured, or if the system is calibrated with a 100% reflector 


in the shock-tube, the true reflection coefficient of the shock- 
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front may be determined. 
5./.2 Spatial Resolution 


If the output signal shows no signs of beat amplitude 
degradation due to curvature of the shock front, values of the 
reflection coefficient across the entire horizontal diameter 
of the shock front can be calculated from each cycle of the 
waveform giving a lateral spatial resolution of 2-3 mm. The 
measurement at any single point is made in a time during which 
the shock wave moves about 1 mm down the tube. The resolution 
of the measurement in the direction perpendicular to the plane of 
the interferometer and shock tube is at the worst 1 mm, and 
considerably less for points near the axis of the shock tube 
where the two beams are brought to a common focus. If the shock 
front is curved reliable reflection coefficient measurements 
can only be obtained for points on or near (within 2-3 mm) of 
the shock tube axis, but the same degree of spatial resolution 
in all three dimensions is obtained. The measurement time is 
in the region of ne sec, and successive segments of the shock 
front are then measured at intervals of rear sec, the whole front 
being scanned in perives sec. It is unlikely that the reflection 
coefficient varies much in this order of time interval, so the 
observed waveform can be used to determine a valid spatial 


variation of reflectivity across the front. 
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5.8 Effect of Shock Front Tilt. 
eee ee POUCA Cs 
5.8.1 Effect on Shock Velocity 


If the measurement is performed at large angles of 
incidence, eq.(5.10) applies rather than eq.(5.9);, and the beat 
frequency is independent of the actual angle of incidence. 

The effect of tilt of the shock front would be to alter the 
effective angle of incidence, so provided the tilt is not so large 
that the reflected beams are deflected entirely out of the aperture 


of the output lens, there is no effect on the velocity determination. 
5.8.2 Effect on Shock Curvature 


Eq. (5.12) shows that the value of the radius of curvature 
varies roughly as a and is thus strongly influenced by any tilt 
in the shock front, since this produces a variation in Oe. However, 
the tolerance of the system to shock tilt can be reduced to such a 
point that only shocks with tilt angles, $¢, much less than on can be 
observed. For a system where the aperture of the output lens is small 
compared to its focal length, F, (this is generally true, since the 
maximum aperture required is just equal to the interaction length 


for a planar shock front, and this length,1 , is given by | = 2x, Cor 
Misa) 


and is thus small for large angles of incidence, i even if apertures 
several times this minimum value are used, to accomodate curved shocks) 


it can be shown that the maximum tilt which can be accomodated is 


Io = cet 36 69 =} gnonsypansm ‘et? ae 
Oe, 
{© 2) .08 Nels youtss Asti OTR. pe. 


se a 

haat yofehs 'audon Sid io shobhegenee e 
eit 

ait soste pi. od bluow gnet: dqiore std To 4it9 ai a 


, Ee: mia 3! 10 603 omiotiiomt 20 sigan” 


hoe 
- 


duo Vist! ‘sh ors empsd tosustie® oil: 
af | 
ze ' : hy v4 « \* ; i ‘ as Of me | awe | orf , erat uqau0 ‘i a 


’ ba ty 


rr GeV I Houtte pe : 326: a% 


— ee a 


oy wine Fu ia ett } } 30 ¢3 Ant? avon fe, (SUR oe we 


. ; ated. adie nisi 
‘4 =? ~ } 
» 
/ 
& 09 ) Beas i , ,50ne 3°10 JAN = “ng =f 4 


i’ ; 
ad eno 6 nal? west! doum ,4 feslgne 2013 dotw exacts vine ant 


: Fv. ; 

i - oy | = . ° . A - 7 

Linwe 2} -easl fuqdeo Silo To sivimsqs 53) ezede er sr aA + 
WwW a 

s@3 eonia oust yilevenca et eidg) 7 “WHA8 asat leant aah oD 


ftanst nor} [793i odj <P Leugo= Jeu rE at bantu es 


I ao 38 = ; vd navigres , s figest aids eae mena 


(€1.2) + _ ves 3 
- 


re] 
eanvivone ii wave }. »sonebronl. jo eo lanes “Opae 


. ni ¢ } . = ; ; 
eras atebomd aie G2 , bear oie su Lay cliabie f 
7 & ania 


ak tatmduemsnn ott. aD totete ahaa 


oe 


- a 
eae. 
7. 


em) 


given by 
ato (5.14) 
ee Pee 
cot i 
or > = 20. ya = Pie (orl) 


So $/6 can be reduced to any value required by using lenses with 
very long focal lengths. 


In the system used, 


Ce = 25 mm 
F = 250 mm 
so o/6 = 0.1 


and a maximum error of 20% in rs can thus be expected. 
5.8.3 Effect on Reflection Coefficient 


The reflection coefficient is expected to vary roughly 
as rt (where 8, = I 1) (eq.3.45), and for the values of 8» 


used here, we have 


-4 -4 | 
(Ge ig) Ban 0 oe(ln Alas) (5.16) 


so that for a value of $/8, of 10 al (section 5.8.2) a maximum error of 
40% in the value of the reflection coefficient would arise if the shock 


front was tilted. This uncertainty is, of course, inherent in any 
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determination of reflection coefficient. 
oS. Experimental Results 


Fig. 5.2 shows the beats observed experimentally with 
the interferometer shown in fig. 5.1. Fig. 5.2.1 shows the beats 
obtained with a plane mirror in the shock tube, with an aperture 
stop placed at the photomultiplier to give a 100% modulation, as 
described in section 5.4. In fact only 85% modulation is achieved, 
due to stray light, and the splitter reflectivity differing from the 
optimum 38% (36% or 41% reflectivity would reduce the modulation to 


85%). 


Bigs 5.2.2 shows the beats observed with a shock wave, all 
three of the major recombined beams being permitted to enter the 
photodetector, to give the widest possible system aperture. Under 
these conditions a maximum modulation of 60 to 70% is expected 
for the range of reflectivities which give 85 - 100% modulation as in 
fig. 5.2.1. In fact 50% modulation is observed, over the central 
part of the waveform, indicating the areeeee of some stray light, 
Orewa aeit jd eteemene in the optical system. 


The relevant parameters for fig. 5.2.2 are:- 


Vertical Sensitivity: 2.6.x ORE w/cm (both beams) 
zero lines displaced 2 cm. 
Time Scales: 0.2 usec/cm and 10 usec/cm. 
Shock Velocity, v: (46844)m/sec [Mach 1.4) in air at 720 Torr and 


293°K. 
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FIG. 5:2 BEATS OBSERVED WITH DOPPLER INTERFEROMETER 


FiG.521 BEATS OBSERVED WITH PLANE MIRROR IN TUBE 
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Angle of Incidence, i, = 87.82° (40.02°); 0=2.78°. 


Angle between Beams, 6i = 0.58° (+0.02°). 


There is Seri indication of shock curvature at one side of the tube 
since the beats die away towards one end of the pulse, but no value 
of the curvature could be determined in this case, since beats were 
present right up to the other edge of the pulse. The minimum value 


of R would be 300 mm, 
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FIG.5:4 BEAT FREQUENCY VS. BEAM SEPARATION 
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The results of a systematic investigation of the validity of 

eq. (5.10) are summarized in figs. 5.3 and 5.4. Each point on the 
graph represents the mean of at least five determinations of each 
quantity, es identical experimental conditions, and the error- 
bars represent the experimental standard deviation. No variation 
in frequency with angle of incidence was found, as predicted by 
eq. 5.10, at the large angles used, 80-88") . An output pulse 
was observed, but with no beats present, if either one of the 
incident beams was obstructed. 

The shock velocity was varied by using different diaphragm 
thicknesses and driver gases, and by varying the downstream pressure. 
It was found that reliable velocity measurements could be made for 
downstream pressures down to about 80 Torr, below which there was 
insufficient reflected light to give accurate measurements. At the 
lower pressures significant values of shock curvature became 
noticable. Fig 5.5. shows the type of effect observed, the two traces 
shown being calibration traces taken with a plane mirror and a 140 mm 
radius curved reflector. Individual beats are not resolved in 
these measurements, since a large value of 8 was used, to give a 
reasonable total pulse length, with the relatively small radius of 
the curved reflector. Fig 5.5.3 shows the effect with a shock wave 
as the curved reflector, giving a radius of curvature of (6.0+0.4)m, 
for a shock wave moving into air at 720 Torr at a velocity of 
(635+5)m/sec [Mach 1.90]. Fig 5.6 is a summary of five such values 


determined at each of several downstream pressures, Py? showing also 
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IN SHOCK TUBE 


PLANE MIRROR 


FIG.5-5 SIGNALS OBSERVED WITH CURVED REFLECTORS 
FIG. 5-5-1 
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FIG. 5:6 SHOCK FRONT RADIUS VS. PRESSURE 


FOR AIR AT MACH NUMBERS LESS THAN 3 
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the theoretical dependence of Ron Py predicted by Duff & Young SO 
and a representative value taken from the results of Linzer & 


La 
( e Mazzella & de Boor eae working in Argon have obtained 


Hornig 
similar = Re at 1 atm but with considerably lees accuracy. 

Their results differ by a factor of 5 or so from those reported here 
at pressures less than 1 atm. It is possible that these discrepancies 
between theory and experiment are due to the surface roughness of 


the inside of the shock tube (11) 


which was machine finished to 
~ axl Or mm, for the results plotted in fig 5.6. 

Fig 5.7 shows the experimental variation of (coherent) 
intensity reflection coefficient with angle, together with the 
theoretical dependence predicted by eqs. 3.33 and 3.41. Good 
agreement is obtained, over the full range of angles used (a to 6° 
grazing incidence). The experimental points were obtained with a 
Mach 1.90 shock wave, using an hy tate et angle é6i (or 66) of FS 
Each point on the graph represents 5 or more measurements, so that 
the effect of random variation in shock tilt which might be expected 
to produce a 40% error in R is reduced, giving an error of 10-202 
for each point. A comparison of fig 5./ with fig 4.5 shows 
immediately the improvement in agreement with theory given by the 
measurements taken with the differential interferometer, particularly 
at 57 <6F s It is evident that there is normally a large amount of 


scattered light present at these angles - as much as ten times the 


true reflected signal, as given by the interferometer. 
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FIG.5.7 REELECTION COEFFICIENT VS. ANGLE (EXPERIMENT) 
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VI Alternative Measurement of Shock Curvature and Tilt 


6.1 Introduction 


In chapter V it was shown how the curvature and reflectivity 
of a shock front could be measured to Bmract es of 20% and 40% 
respectively, the errors being due to the possibilities of tilt in 
the shock front. In this chapter an alternative method of determining 
curvature and reflectivity, which also detects shock tilt will be 
described. This experiment also is capable of detecting the presence 
of long-wavelength (i.e. ~ shock tube diameter) perturbations in 
the shock front of which shock front curvature is merely a special 


(Ist order) case. 


6.2 Experimental Arrangement 


The arrangement of the optical system for these measurements 
is shown in Fig 6.1. The incident laser beam is focused down at a 
single interaction point on the shock tube axis, and light reflected 
from this print by a plane shock wave, perpendicular to the tube axis 
(i.e. zero tilt) is obstructed by a 2 mm opaque strip, mounted across 
the photomultiplier aperture. Angles of incidence between 70° and 38° 
are used, as before. 

As a plane shock front passes through the test section, a 
square-pulse is observed at the photomultiplier output, with a dip 


to zero near the pulse CenvceL,; corresponding to the obstructed 


axial reflection. The (constant) pulse height gives the reflection 
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coefficient of the shock front, as in the preliminary experiments 
(chapter IV) and the position of the zero-dip indicates whether the 
shock front is tilted, relative to the shock tube axis since the dip 
is alg exactly central for reflection from the mid-point of an 


untilted shock front. 


Any variation in pulse height during the reflection process 
esethen due to the reflection Occurrinpat different angles of incidence 
as the shock front sweeps through the interaction region - i.e. 
curvature of the shock front - giving a steady increase or decrease 
in reflected aeeaee across the pulse for concave and convex shocks 
respectively; or a long-wavelength perturbation of the shock front 


giving a more complex variation in the output pulse height. 


The analysis of the output signal from a generalized 
shock perturbation is extremely complex, but for the most interesting 
and experimentally more common occurrenceof a simple shock curvature, 
the radius of curvature can be calculated relatively easily from 
the ratio of the maximum to the minimum pulse heights. The 
possibility of shock tilt invalidating the results is eliminated 
here, since measurements are only taken from those output traces 
with a central ‘dip’, indicating zero shock tilt. The true 
reflectivity of such shock fronts is determined by extrapolating the 
trace into the central 'dip' region to obtain the reflected intensity 


for the on-axis point, for which the angle of incidence is well-known. 
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6.3 Determination of Curvature 


Consider a shock front having a radius of curvature, 
to: propagating down a shock tube of radius ree Then a light beam 
reflected off the shock front at a nominal grazing angle 8, will 
in fact be reflected at an angle which varies from (8 +0) to 


(6 -a) (fig. 6.1) where 


es 
7 t 
og =m (6.1) 
Ss 
The intensity reflection coefficient, R, is given by 
ede (5.45). viz. 
1 6 = it on 
R= Tate 1 Salles) (6.2) 
46 6 
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FIG. 6:1 CURVATURE MEASURMENT 
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8S = 6-a) to the minimum (6) = 6+0) is:- 
1 one qt 6 
I 4 Ce) [1 2 e) 
max _ 4(0-a) (6-a) 
Tay 
abe 1 én. 2/ 1 én 
eR. Oe rs ry. 
4 (6+0) (8+a) 
1 6n 
C=) 
" sin p= oa)" 
-O 
l1-—+,¢4 
Z n 
(8+a) 
max 0+ 4 
or rit = Bea | to the first order 
min 
qo 
so that j/_max 4 
ANY an 
a= 86 ae 
at tbr 
NCR 
' “min 
sis 
r (6.1) 
s 
oe ite 
giving _ max, 4 
r 4 ; 
Py OE 
ts wind (jon, 
ik 
~\ /_ max l 
Ai 
min 


The above derivation ignores the effects of the finite size of the 


(6.3) 
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light beam, and the exact beam profile. These effects have been 


205 who derive 


dealt with in a more exact analysis by Pert & Smy 
an expression which reduces to eq.(6.7) provided that the effective 


light beam width, d, is small:- 
seg te = ea (6.8) 
: ; 


For the parameters used in these experiments this yields the condition 
that 


TA 500 mm (6.9) 
a condition which is found to be satisfied in practice. 
6.4 Determination of Shock Tilt 


If the shock front is tilted through an angle, 9$, the 
dip in the trace due to the small obstruction at the detector will 
not occur at the mid-point of the reflected pulse, as set up in 
the initial alignment procedure, but will be earlier or later, 
according to the direction of the tilt. 

The geometry of the system is shown in fig 6.2. The dip 
in the signal occurs when the shock is at BE, rather than at AE 
i.e. when the shock is an additional distance, X, downstream of the 


centre (alignment) position. 
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G62 TUT MEASURMENT 


TILTED SHOCK 
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clearly, x = ue ét 


(6.10) 


where 6t is the time difference between the centre of the reflected 


pulse, and the observed dip in the signal 


also AB = —— 
sind 
D 
oa cosé 
so in A BAC 
peg xn Uoper Limit Codie ras 
sin® sin 26 ~~ cos®@ sin 2(0-6) 
or x = D tanég goer? 


sin 2(6-6) 


so for small angles 6: $9, 


2 | > 
x ~ Dé Ce 
\ 
x0 
we re 


(6.11) 


(6.12) 


(6.13) 


(6.14) 


(6.15) 


(6.16) 


(Ose) 


The maximum value of x for which reflection can still occur 


is limited by the shock tube radius, and the angle 6, to 


vr a ep | 


ar. id 
‘ : 7 4 i] 
bd os ae 

: > - 

ly 8 
| rie 
, st oh 7 vit 
6 ea 


it! 


7 


. Pre 
‘ Sosngzattih ants ada rs 


cit, hevaobey eff 


Ne 
T 
. aa 
a] 
os 
NN 
& 4 
: 4 
G | An 
\ 
. parr 
| ~{OL.A) po vd a 
e 
a6 Y¥ 
(Vi. >) 3 
36 yw + Of 
’ ; 
: ; 7 
swuase [123e mao oltrelio: doin 297 x to sviev mae 
7 ay ) : 


o3 <f eighe ois bra ,2ulbber o¢ 3 Anod 


ee 


78 


max iE Ss max 


so we have 


o = 
max 
Dot a 
and normally Bae. 
rt 
So. = —— 
4 nese ; D 


giving an upper limit to the measurable tilt of about 5 m.rad 4°) 


(6.18) 


(6.19) 


(6.20) 


aha 21) 


in these experiments. In fact the observed shock tilt was always very 


much less than this, typical values being + 0.5 m.rad. 
Under these conditions (6.17) and (6.20) can be used to 


derive the more convenient approximation:- 


which is sufficiently accurate for most of the experimental 


situations encountered here. 
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VII The High Intensity Shock System 
7.1 Purpose 


In chapter 3 it was shown that the practical upper limit 
to the Mach number produced in a conventional pressure-driven 
shock tube is about 3.5, and in fact it is inconvenient to work 
much above mach 2. A shock tube was therefore required onien 
would produce the much higher Mach numbers - in the region of 


(24) 


Mach 10 - which are more often encountered in plasma - dynamic 


situations. 
7.2 Production of Strong Shocks at High Pressures 


wright ‘1°? hasreviewed the methods by which strong shocks 
can be produced. With the exception of explosively driven shocks, 
however, it has only been possible to produce strong shock waves in 
gases at low pressures (i.e. a few Torr), the most common form of 


(25) 


apparatus being the electro-magnetically driven shock tube or 


the combustion shock pape ts02 


More recently, electrically-driven 
shock tubes have been designed which combine the ease of operation of 
the electromagnetic type, with the absence of circulating currents in 


3 
the plasma created behind the shock ppc a 


In this type of shock tube, a high energy capacitor bank is 
discharged through the gas in the driver section of a diaphragm-type 


shock tube, thereby producing intense heating of the driver gas, 
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giving rise to an extremely high sound speed, and a very high pressure 
ratio across the diaphragm - both conditions being essential to the 


(19) 


production of strong shock waves (section 3.1.7). Phillips has 
shown that the high Mach numbers produced in this type of shock tube 
are due only to the intense heating produced by the discharge, the 
electromagnetic forces involved in the production of the shock wave 
being negligible. The term "electro-thermal" is therefore used to 
describe this type of shock tube. Its use has so far been restricted 
to producing shock waves in low-pressure gases. There is no reason, 
however, why extremely strong shocks should not be produced in gases 
at high pressures (a few atmospheres) by an electrothermal shock 


tube which is designed to operate at the extremely high temperatures 


and pressures which are required, thus extending the range of previous 


work 619» 30) | 


7.3 Theory of the Electrothermal Shock Tube 


From Chapter 3, (eq.3.25) the diaphragm pressure ratio can 
be numerically related to the shock strength, y, for a Helium-Air 


shock tube by: 


a y ein) 


C. P : 
TO TG 20g eee 


where C_ and C. are the sound speeds in the shocked and driver gases. 
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The Mach number is given by (eq. 3.17):- 


a ea 


- ‘by +1 
M So (jez) 


from Eq. 7.1 it is apparent that even with an infinite diaphragm 


pressure ratio we still have a finite shock strength, given by:- 


C 


0.629 — 232 2 i (3) 
3 Voy +1 
| Cy 2 
or en eee Lyon. (7.4) 
3 
giving AS ley jBPae, 
and M = 5.8 
max 


for Helium-Air at 293°K. 


If we assume ideal gas behaviour, then if the driver gas temperature, 


T35 is raised from room temperature T3099 we have 


Cy SNe ry T, Co) 


oO rf a 
30 30 
p aL 
and hy = (7.7) 
Pa0m nO 
where C43 P39 are the room temperature sound speed and pressure of 


the driver gas, and ty is the gas constant for unit mass of the driver 
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gas. 


Equation 7.1 then gives:- 


eS 6 e ek eee ee dae (7.8) 


for Helium-—Air. 


Fig. 7.1 shows this dependence of driver gas temperature (and 
pressure - eq. 7.7) on shock strength and mach number for a system 
with an initial filling pressure P39 of 5 atm. (Helium) and air 
at one atmosphere in the downstream section of the shock tube, the 


entire system being at 29351 initially. 


From fig 7.1 it is apparent that although the gas temperatures 
required to obtain Mach 10 shocks are extremely high, they are 
within the range which can be produced by intense discharges. The 
pressures involved are also considerable but are well within manageable 
limits. It should be realised that eq. 7.8 and fig. 7.1 represent an 
upper limit to the Mach numbers and shock strengths which can be 
obtained, since no allowance has been made for the variation of Y3 
with temperature (as the driver gas becomes ionized)or for any 


departures from ideal gas behaviour at high pressures. 
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FIG. 7:1 DRIVING CHAMBER TEMPERATURE AND PRESSURE _VS. 


MACH NUMBER AND SHOCK STRENGTH - 


FOR He-AIR, INITIALLY AT 5 ATN.,1 ATM. 
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The actual Mach number can be predicted for any given shock 
tube if the enthalpy of the driver gas is known. If we assume that 
the entire energy input W (stored in the capacitor bank) is 


transformed into thermal energy of the driver gas we have:- 
BATS 
P, =¥ (7.9) 


where V. is the volume of the driver chamber. 


3 


In fact there may be considerable energy losses in heating the 
walls of the chamber, and in vaporizing material from the electrodes 


and diaphragm - more exactly 

Per eg 

= ta (7-10) 
where ¢« is defined as the "thermal efficiency" of the shock tube. 


Phillips "2? found values of e« of about 0.45, in his shock 
tube. Assuming a generally similar configuration we are now in a 
position to determine the energy input per unit volume of driver 
chamber required to produce a Mach 10 shock wave, as about 100 Joules 
for every cubic centimetre or energies of around 10 KJ for a volume 
of 100 cc at 5 atm. These figures are quite easily obtained in the 
laboratory, and the construction of a high-intensity shock system 


to operate at 1 atm, downstream is thus quite feasible. 


7.4 Design of the Shock System 


An existing capacitor bank was used to provide the required 
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energy. This bank had a capacitance of 43 uF at 20 kV, giving a 
stored energy of 8.7 KJ. The volume of the driver section was 

fixed at 250 cc, giving a maximum driver pressure at 100% thermal 
efficiency of 340 Atmospheres, and a maximum mach number of 11.4. 

The driver section was designed to withstand a maximum internal 
pressure of 680 Atm. (10,000 p.s.i), and to give an additional safety 


margin over the normal pressure-vessel design figure. 


The downstream part of the shock tube must be capable of 
withstanding the shock overpressure, given by the shock strength y, of 
150 Atm. This is well within the limits of normal high-pressure 
Steel tubing, but in order to accomodate wall-mounted piezo-electric 
pressure probes and ionization probes, a much heavier section steel 
tube was used - 32 mm o.d. by 15 mmi.d. This gives sufficient 


length to a standard fine-pitch threaded plug to withstand 150 atm. 


A longitudinal section of the shock tube is shown in fig. 
7.2, stainless steel being used throughout the system, with the 
exception of the insulating backplate, which was made of “rerloni. 
The driver section was made conical in form to obtain the maximum 
flow velocity of the gas from the driver, after the diaphragm had 
opened. This has been shown to give a 10% increase in Mach Ae 
and also assists in maintaining the shock for as long a time as 
possible. The driver chamber length of 100 mm was sufficient to 


give a reasonable (170 mm) useable length of shock tube. For 


distances greater than about 1.7 driver chamber lengths (under the 
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FIG. 7-2 HIGH INTENSITY SHOCK TUBE 
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conditions present after the discharge) down the shock tube the 
rarefaction wave reflected off the back of the driver section ol 


"catches up' and interferes with the shock lieve bode 


The diaphragm is retained between two 25 mm diameter O-ring 
seals by means of a quarter-turn single start interrupted screw 
thread, so that diaphragms can be changed quickly and conveniently. 

A number of different diaphragm materials were tried. The most 
important properties required are that no fragments of the ruptured 
diaphragm should be driven down the shock tube to perturb the shock 
wave and damage the instrumentation, and that the diaphragm should 
not rupture immediately the discharge is initiated, but after most 

of the capacitor bank energy has been transferred to the driver gas. 
Two materials were found which satisfied both requirements - 0.025 mm 
thick "mylar" film, which was vaporized by the discharge in a time 
comparable to the total duration of the (underdamped) discharge, and 
0.25 mm thick steel shim stock, with two 25 mm long slots 0.24 mm 
deep milled across it in the form of a cross, so that the diaphragm 
folded back against the shock tube inlet, after completion of the 
capacitor bank discharge. Because of the greater convenience, "mylar" 
film was used for most of the experiments, although the velocities 


obtained were less predictable than with the steel diaphragms. 


The vacuum, purging and backfilling systems for the shock 
tube and driver sections were equipped with electromagnetic values, 


so that the operation of the shock system could be made largely 
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automatic. The whole system could be cycled from shot to shot in 
less than a minute. The discharge was initiated by an additional 
trigger electrode, in the driver chamber, thus eliminating the 


y 9) 


energy losses (up to 50Z inherent in using an external spark- 


gap switch. 
7.5 Operation and Performance of the Shock Tube 


In a Pee cial shock tube, the finite opening time of the 
diaphragm, with rupture first occurringon the shock tube axis, 
results in an initial radial velocity and axial acceleration of the 
shock front, the maximum axial speed and optimum stability normally 
being attained at a distance of 10 to 20 shock tube diameters 
downstream of the diaphragm. For the shock tube described here, 
using "mylar" diaphragms which were vaporized completely by the 
driving discharge, this effect was considerably reduced, and it was 
found possible to work within a few diameters of the diaphragm. Most 
of the measurements were taken at points from 8 to 10 diameters 


downstream, however. 


Piezo-electric pressure transducers mounted in the shock 
tube wall were used initially to measure the shock velocity, but the 
ceramic piezo-electric crystals were rapidly destroyed, apparently 
by the sudden reduction in overpressure as the rarefaction wave 

(i5) 


followed the shock wave down the tube . Sufficient results were 


obtained, however, to compare with velocities obtained with ionization 
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probes, to check that the ionization front immediately behind 

(and caused by) the shock wave had the same velocity as the shock 
front. Modified chassis-mount B.N.C sockets were used as ionization 
probes, screwed directly into the shock tube wall, so that the 

probe was flush with the inner surface of the tube. Three probes were 
used at 25 mm centres. Voltages were applied between the upstream 
probe and each of the downstream probes, so that current pulses 

were generated as the shock (or ionization) front passed the second 
and third probes. A typical oscilloscope trace of the voltages 
generated across suitable load resistors is shown in fig. 7.3. The 
traces are initiated by the capacitor bank trigger pulse, and show 
the discharge current in the form of "pickup". By extrapolation of 
the shock velocity it can be seen that the diaphragm opens when about 
90% of the capacitor bank energy has been transferred to the driver 


gas. 


Fig. 7.4 shows the measured Mach numbers produced by the 
shock tube as a function of capacitor bank voltage, with the theoretical 
lines, for a number of different values of « (eq.7.10). It was found 
that ¢ decreased after a number of shots and this effect was traced 
to the gradual contamination of the insulator surface by evaporated 
metal and diaphragm material, resulting in 'tracking' and giving rise 
to a less efficient form of discharge. The minimum‘value of ¢ was 0.3 
which is reasonable, and values very close to one were obtained 


immediately after cleaning the insulator. A large number of different 
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FIG. 7-4 CAPACITOR BANK VOLTAGE vs, MACH NUMBER 
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insulator materials were tried, but no other material was found 
which had such a good tracking resistance as "teflon" coupled with 


the required mechanical strength. 


The performance of the shock tube was generally similar 


to that of the much larger system of Camm & Rosa Ole which was 


operated with downstream pressures of 0.1 to 0.3 Torr, and somewhat 


better than that of Bh Misty suomi 
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et Measurement of Reflectivity and Curvature of Hypersonic Shocks 
8.1 Reflection Coefficient 
8.1.1 Choice of Method of Measurement 


A number of additional experimental difficulties are 
introduced in the measurement of reflection coefficients of shock 
waves having very high Mach numbers. At Mach numbers greater than 
about 4, the gas behind the shock front is heated sufficiently by the 
passage of the shock wave to cause the gas molecules to become eer edaue. 
and to emit light in the visible region of the spectrum. Ultimately 
the gas becomes ionized, at very high Mach numbers. The light emitted 
by the shock heated gas is intense enough to mask the reflected laser 
beam, even at extreme grazing incidence, and the situation is considerably 
worse at larger grazing angles with the shock tube used here, since the 


photodetector system can also see the plasma produced by the driving 


discharge, as it propagates down the shock tube behind the shock front. 


The presence of this large amount of 'stray' light immediately 
suggests the use of the Differential Doppler Interferometer described 
in chapter 5, and used successfully to measure the reflection coefficient 
of low Mach number shocks (section 5.7). However, the frequencies 
generated by shocks moving at Mach 10 or so would be in the region of 
100 MHz or more, for grazing angles up to about 10 degrees, and would 
thus be outside the maximum bandwidth of any state-of-the-art detector 


system for the expected reflected intensities (chapter 2). At angles 
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very near to grazing incidence, the frequencies and reflected intensities 
produced might just be manageable, but only about 10 fringes would be 
observed, even in the complete absence of any curvature or tilt of the 
shock front, due to the extremely short interaction time (0.1 usec) 
between the shock wave and the laser beams which in turn results from 

the high Mach numbers involved and the small diameter (15 mm) of the 
shock tube. In practice, of course, a shock tilt of around 2 degrees 

or the equivalent curvature (400 mm shock front radius) could very 
easily be present, and it is unlikely that more than two or three 

fringes would be observed - insufficient to give a reliable reflection 


coefficient measurement. 


It was therefore decided to use the simple method of measuring 
reflection coefficient described in chapter 4 (section hd ls) eS 
much stray light as possible was excluded from the detector system by 
carefully screening the entire optical system and interaction region 
from the room lights, and by using an interference filter centered at 
632.8 nm wavelength, to exclude as much as possible of the light 
produced in the shock tube. Filters with bandwidths of 20 and 2 nm were 
used, and it was found that there was little to be gained from using 
the narrower bandwidth since this filter also attenuated the laser beam 
considerably (x0.3), and thus only slightly increased the "signal to 
noise ratio" of the detector system. The normal Doppler shift of the 
reflected laser beam is of the order of Ageeian so there is no 


possibility of the wavelength of the reflected beam falling outside the 
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pandwidtheof-the filter. 
8.1.2 Experimental Results 


The measurements were made as described in chapter 4, at 
points 75 to 150 mm downstream of the shock tube diaphragm. The 
photomultiplier was operated with a 50 % load resistance to give 
sufficient system risetime to accurately observe the reflected pulse 
whilst retaining adequate output voltage, at the grazing angles of 
incidence between 4 and 10 degrees which were used. Angles less than 
4 degrees were not used since the laser beam had to be stopped down to 
an impracticably low diameter, to avoid reflection of only a part of 
the beam occurring. Above 10 degrees, the background light intensity 
became too high and the reflection coefficient too low to observe any 
meaningful reflected pulse, even with a 2 nm notch filter ab the: phote— 
‘detector. Fig. 7.3 shows a typical measurement, the upper trace being 
the reflected light intensity, and the lower trace the ionization-probe 
signal for velocity measurement. Fig 7.3 Leowtorea machi .2 Suock, 


giving a reflection coefficient of faders LO 


Fig 8.1 summarizes the results obtained at a number of 
different angles over the full range of Mach numbers available. The 
results obtained at velocities below mach 2 or so are somewhat 
unreliable, due to the low reflected light levels observed, comparison 
with figs. 4.5 and 5.7 showing that the results obtained in this region 


are as much as a factor of 3 different from the earlier results. The 
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FIG. 8:1 REFLECTION COEFFICIENT VS. MACH NUMBER 


FOR AIR AT 1 ATM. 


10 210 310 (m/sec) 


: 
j< . a “ay ‘Tif ss 
‘i + 8 A ri By x 
e > @ eee w= - - 7 La 


‘ 
/ 
( 
_ - 
‘4 ¢ ; 
‘ 
- 
-; 
4 
1 


a 

Tat ae . Le a md me 

Ss ol 7 “Sener, 
; ; + nN 


7 


discrepancies are random however, and the general trend of the 

results is in reasonable agreement with those obtained earlier, and 
with theory. Fig. 8.1 also shows the theoretical variation of 
Reflection Coefficient, Ry with Mach number at a number of different 
angles of incidence, appropriate to the experimental points, the 
curves being calculated from eqs. 3.33, 3.55 and 3.56. It will be 
seen that the experimental values of reflection coefficient are in 
good agreement with theory (Ry = R, at these angles) over the entire 
range of angles and Mach numbers, with a possible trend to higher than 
theortetical values, at the lower Mach numbers (2-4), at all angles. 
It would therefore seen reasonable to assume that the shock reflection 
coefficient follows the theoretical variation with angle and Mach 


number up to the values associated with, for example, laser-produced 


plasmas (M “30, 8 45°), 
8.2 Measurement of Shock Curvature 


The curvature of the shock fronts produced by the high-intensity 
shock system were measured by the method described in chapter 6 (section 
6.3), except that the central Jae ee strip used to detect shock 
front tilt was not used, because of the short length of the reflected 
pulse. In view of the large but consistent curvatures (measured at 
various angles of incidence) it is apparent that any tilt in the shock 


front is randomly oriented, and of under 1 degree. 


The values of the radius of curvature of the shock front were 


determined from the same traces as the reflection coefficients, and are 


- . 
a ee 


} Yo keess IetedSg Si) bin , Tavswan “ob ey 


G Ls 
LO fin tansy : sroSsr 3s ent i i [s [3 -ait 
a 
imeyetI£b to x un & dn soduur dost dsitw A ‘giobobiiedD lan 


. : inn i} 4 
oa3 pintea Latasmii19qgxs of } Mik o7q7g6 a Mohs ont 9 as. 
” “ - 
4 


_ Z " # 
iw ot .OR,€ bab 28.E .f2,E .ap } bevnluchso giited: cov 


Jnolotiises nokgastist to « 1 LetaoadaoGes 982 ‘ses oa 
iw 

Sa 7 “ — “ a “ 
“ee & 9 im 36. , 1) yraots AIIw 1omaa ge boog 
j © 


hd of baans ofélesoq 5 div .2ersdapin fost bas eaiges to op 


; : . ' ep yer va 

eolene Iie ta (A-)) avedawa « ws $07.38 ,3eauiey fxoties 198d? 
; ? ioe 

ipeltey dooms 903 fnilt ouveze 03 ofderiesey qpoe SraleeMe ea 
ise! ban olgnc djlw noltetwy Joolsavoads-SNg nwo Litto’ snot sty ie¢ 


besuhotq-veest ,siqwexs ro? ,130 ig.tcoee oulsvy od) 62 @y 1° 


. g o 
{90/2 %¢ Igemaxnugest” § 


vyiteqsenl-igit sfd <d bssvbor, ox? A502 arid lo sayvésavebe sat 
/ : a 
i 
eobyose) 8 tovantio at baditseab bolton afd yd bsavebem ayew novege json 
a + % 


4 j bd f. we 4 ‘4 

dsede zosteb of bsauv qtrte aoaitouriado. I[essneo sid 3683 fapae a 

beroottex a2 Yo digaal atode af3 to setgosd ,bsett. an) eow why 1107 
| ‘ 3 7 ; » - i 

29 botuaesm) eoiwisvips Jpase leno tnd agubl sia to wa. al 98, 


é TS : 
dnote odd nit 31%3 yas seid Inetagys etouit (ssrsbloah te a 


5 ana 
F : wer 


Ea “pered { reheat to bos. gbesnetie- 
ie : 7 4 ia a) 


“row $082 docs ait i. sxvievaus. io gests i sis. 2g > Bol 


eit 

ants 2?909 mokrowtter 
., i 
; 7, Fa | 


98 


FIG. 8-2 SHOCK RADIUS OF CURVATURE VS. MACH NUMBER 
FOR AIR AT 1 ATM., 150mm. DOWNSTREAM 
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shown as a function of Mach number in fig. 8.2. It is apparent that 
there is no correlation between the shock radius and the grazing angle 
of incidence, as expected, and the results are replotted without 
distinguishing the JET Oteincidence inelie. 6.35, which indicates 
that the shock curvature (= ~) is directly proportional to the Mach 

Ss 
number, or shock velocity. It is noteworthy that there is considerable 
"scatter' and in the results at low Mach numbers (<5) where all the 
previous work on shock waves at atmospheric pressure has been carried 


out, so that the linear increase of curvature with Mach number can only 


be detected using the High Intensity Shock System. 


Fig. 8.3 also shows some curvature measurements taken at a 
different point along the length of the shock tube, indicating an 
approximately linear increase of shock front curvature with distance of 
propagation down the tube. This behaviour is consistent with the growth 
of a boundary-layer at the shock tube walls, as the shock wave 
progresses, the boundary layer in turn causing disturbances which 
distort the shock front, to a theoretically parabolic profile, with the 
vertex of the parabola touching the shock tube eeu inp practice, 
the measurements described here are too approximate to determine the 
exact profile of the shock front, and it is more useful to invoke the 
concept of the maximum extent, or maximum protrusion, x of the shock 


from a plane surface. For aspherical shock front, it is easily shown 


that for small values of x 
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where tS is the shock front radius 


and r, is the radius of the shock tube 


x is thus an alternative measure of the shock curvature 


(l/r) and is plotted as such in figure 8.3. 


(28) 


De Boer, in an involved but still approximate analysis 
has depeunned the theoretical variation of x with Mach numbers for 
both laminar and turbulent boundary layers, in a shock tube of 
circular cross-section, the turbulent boundary layer being appropriate 
to downstream pressures greater than 100 Torr. These theoretical 
curves are shown in fig 8.3, which shows quite clearly that the trend 
indicated by the curve for the turbulent boundary layer is followed 
by the experimental results, as expected, although curvatures 
considerably greater (i.e. smaller shock front radii) than predicted 
were observed. Indeed, the experimentally observed curvature appears 
to be still increasing, even 150 mm downstream of the shock tube 
diaphragm, and it appears unlikely that the boundary-layer (and hence 
shock curvature) ever grows to its equilibrium value (which seems 
to be at least a factor of three greater than that predicted 
theoretically), since the shock rarefaction catches up with the shock 


wave some 170 to 200 mm downstream of the diaphragm. 


It is possible that the enhanced curvature observed 150 mm 
downstream is due to the proximity of the rarefaction wave affecting 
the boundary layer behind the shock front, but since the contact surface 


is observed by laser reflection at some angles of incidence, to be 
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about 20 mm behind the shock at this ee it would seem unlikely. 
It is more probable that the theory becomes very approximate for such 
a small-diameter shock tube, and a full investigation of shock 
curvature effects at igh pressures must await the construction of a 
larger diameter shock tube, capable of operating in the same régime. 
It may also be remarked that the turbulent boundary layer did not 
produce any marked shock attenuation - the shock velocity remaining 
appreciably constant (to within 10%) over the first 150 mm of the 
tube, under any given set of beat tonke 

The limitation due to finite laser-beam width, given by 
eq. 6.8 (p. 75) yields for this shock tube (radius = 7.5 mm) the 
condition :- | 

t >> 80 mm 
which is easily satisfied for all the results obtained (minimum 


value of Cae 500 mm). 
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IX Conclusion and Suggestions for Further Work 


9.1 Conclusion 


The reflection of a light beam by a shock wave has been 
investigated both theoretically and experimentally, and it has been 
shown that the results observed can be adequately explained by the 
exact application of the laws of optics. Experimental confirmation 
of the validity of the theory has been obtained at Mach numbers and 


ambient pressures which have hitherto been unattainable. 


It has also been shown that the presence of light 
reflected by shock waves is likely to predominate over the required 
signal in a number of plasma-dynamic experiments in which light— 
scattering by the plasma is used as a diagnostic tool, and a method 
by which such scattered and reflected components may be distinguished 
and separated at the light ian has been demonstrated. The 
advantages and limitations of this technique have been investigated 


in some detail. 


Finally, a shock tube ae been developed which is capable 
of producing hypersonic (>Mach 10) shock waves in air at one 
atmosphere, and some new information on the curvature BE such shocks, 
resulting from a turbulent boundary layer at the shock tube wall, has 


been obtained, thereby extending the results of previous workers. 


9.2 Suggestions for Further Work 
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The Differential Doppler Interferometer (chapter 5) 
lends itself to a variety of Plasma Physics Experiments, possibly 
the most interesting of which would be the light reflection and 
scattering occuring at the plasma produced by a eneneed high= 
power (200 MW) laser beam, where very small 'scattering' coefficients 
have been reported, in the presence of extremely strong, (and there- 


fore relatively highly reflecting) shock waves. 


The investigation of these extremely strong shock waves 
at high pressures {one or more atmospheres) would also prove very 
interesting in its own right, if the measurements could be extended 
to grazing angles of 20 to 30 degrees - by the use of a higher power 
laser as a light source - since at these angles the values of the 
reflection coefficient obtained can be used to calculate the gas 
density profile through the shock wave, and information on the gas 


relaxation rates behind the shock front can be obtained. 


The construction of a large-diameter High Intensity 
Electrothermal Shock Tube along the lines of that described in 
chapter 7 would enable more reliable information on the structure 
(and curvature) of high pressure hypersonic shocks to be obtained, by 
the use of a differential interferometer, and would give some 
further information on the effect of tube diameter on shock front 


curvature, 


The present High-Intensity Shock System could fairly easily 
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be adapted to operate with reduced pressures downstream, and with 
gases other than air in the downstream section of the tube, giving 


still higher Mach numbers in a variety of gases, and thus forming a 


(30) 


convenient source of dense highly-ionized plasma It has also 


(29) 


been suggested that a further increase in shock Mach number could 
be obtained by discharging a second energy storage capacitor bank 


into the shock tube, behind the shock front, using the ionization 


produced by the present shock to initiate the discharge. 
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Appendix 


Theory of Reflection from a Beam Splitter 


Consider a plane parallel splitter having a front-surface 
intensity reflection coefficient Ro} and a back surface reflection 
coefficient of one. 

Then a beam of intensity Wy (watt) will be reflected from 
the front surface at intensity 


W, = R 


As Tae (Al) 


and a beam of intensity We will be refracted at the front surface with 


intensity 

hs Ry) Wo 
and will emerge again, after reflection at the back surface with an 
intensity 

WS eR) W, | (A2) 


Then if the initial beams are of equal intensity, 
CATO Wy = Wo PA" mee ee (A3) 


the two main reflected beams will have equal equal intensities 


(as desirable in the differential Dopler interferometer) if:- 


RB, = (1 - R,)” (A4) 
| 
or Re - 3Rp +12=0 | 
giving ae 0.38 (A5) 


as the only root less than one. 


For unequal reflected intensities, Wa? W., the maximum total intensity 


Bp? 


occurs when the two beams are in phase, 1.é. 
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We Shr (Vi, + vit, ) : (A6) 
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A2 
and the minimum intensity occurs when 


W = CoN - i (A7) 


min B 


so the modulation, m is given by 


ae We : Wy bus ee Ae 
a Wy + Ue + 2/W Ww, 


which gives 


(A9) 


and for any observed modulation depth, my the relative values of 

Wa and Wp can be calculated, and hence Ro can be determined from 
equations Al, A2 and A3. It is a relatively easy matter to extend 

the theory to cover the case of multiple reflections in the splitter; 
and for the use of two splitters in succession, as in the Differential 


Doppler Interferometer. In this general case it can be shown ath PEE 


m = 0.85 we have Ry mu).,26.0r 0.41 at each splitter. 
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